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The quasiclassical theory of superconductivity provides the most successful description of diffusive 
heterostructures comprising superconducting elements, namely, the Usadel equations for isotropic 
Green's functions. Since the quasiclassical and isotropic approximations break down close to inter- 
faces, the Usadel equations have to be supplemented with boundary conditions for isotropic Green's 
functions (BCIGF), which are not derivable within the quasiclassical description. For a long time, 
the BCIGF were available only for spin-degenerate tunnel contacts, which posed a serious limita- 
tion on the applicability of the Usadel description to modern structures containing ferromagnetic 
elements. In this article, we close this gap and derive spin-dependent BCIGF for a contact en- 
compassing superconducting and ferromagnetic correlations. This finally justifies several simplified 
versions of the spin-dependent BCIGF, which have been used in the literature so far. In the general 
case, our BCIGF are valid as soon as the quasiclassical isotropic approximation can be performed. 
However, their use require the knowledge of the full scattering matrix of the contact, an informa- 
tion usually not available for realistic interfaces. In the case of a weakly polarized tunnel interface, 
the BCIGF can be expressed in terms of a few parameters, i.e. the tunnel conductance of the 
interface and five conductance-like parameters accounting for the spin-dependence of the interface 
scattering amplitudes. In the case of a contact with a ferromagnetic insulator, it is possible to find 
explicit BCIGF also for stronger polarizations. The BCIGF derived in this article are sufficienly 
general to describe a variety of physical situations and may serve as a basis for modelling realistic 
nanostructures . 



PACS numbers: 73.23.-b, 74.45. +c, 85.75.-d 
I. INTRODUCTION 



The quantum mechanical spin degree of freedom is 
widely exploited to control current transport in electronic 
circuits nowadays. For instance, the readout of magnetic 
hard disks is based on the giant magnetoresistance effect, 
which provides the possibility to tune the conductance of 
e.g. a ferromagnet /normal metal/fcrromagnet (F/N/F) 
trilayer by changing the magnetizations of the two F 
layers from a parallel to an antiparallel configuration 1 . 
However, many functionalities of hybrid circuits en- 
closing ferromagnetic elements remain to be explored. 
Presently, non-collinear spin transport is triggering an 
intense activity, due to spin-current induced magneti- 
zation torques^, which offer new possibilities to build 
non- volatile memories 3 -. Another interesting possibility 
is to include superconducting elements in hybrid cir- 
cuits. When a N layer is connected to a BCS super- 
conductor (S), the singlet electronic correlations charac- 
teristic of S can propagate into N because electrons and 
holes with opposite spins are coupled coherently by An- 
dreev reflections occurring at the S/N interface 4 . This 
so-called "superconducting proximity effect" is among 
other responsible for strong modifications of the density 
of states of N-. In a ferromagnet (F), the ferromag- 
netic exchange field E ex , which breaks the symmetry be- 
tween the two spin bands, is antagonistic to the Bardeen- 



Cooper-Schrieffer(BCS)-type singlet superconducting or- 
der. However, this does not exclude the superconducting 
proximity effect. First, when the magnetization direc- 
tion is uniform in a whole S/F circuit, superconduct- 
ing correlations can occur between electrons and holes 
from opposite spin bands, like in the S/N limit. These 
correlations propagate on a characteristic distance lim- 
ited by the ferromagnetic coherence length -J hD / E ex , 
where D is the diffusion coefficient. Furthermore, E ex 
produces an energy shift between the correlated electrons 
and holes in the opposite spin bands, which leads to spa- 
tial oscillations of the superconducting order parameter 
in F— , as recently observed^—. These oscillations allow to 
build new types of electronic devices, such as Josephson 
junctions with negative critical currents^, which promise 
applications in the field of superconducting circuit o 11 ' 12 . 
Secondly, when the circuit encloses several ferromagnetic 
elements with noncollinear magnetizations, spin preces- 
sion effects allow the existence of superconducting cor- 
relations between equal spins^ 3 -. These correlations are 
expected to propagate in a F on a distance much longer 
than opposite-spin correlations. This property could be 
used e.g. to engineer a magnetically switchable Joseph- 
son junction. These and many more effects have been 
reviewed recentl y 14 ] 15 . 

To model the behavior of superconducting hybrid cir- 
cuits, a proper description of the interfaces between the 
different materials is crucial. This article focuses on 
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the so-called diffusive limit, which is appropriate for 
most nanostructures available nowadays. In this limit, 
a nanostructure can be separated into interfaces (or con- 
tacts) and regions characterized by isotropic Green func- 
tions G, which do not depend on the direction of the 
momentum but conserve a possible dependence on spa- 
tial coordinates. The spatial evolution of the isotropic 
Green functions G is described by Usadel equations 16 . 
One needs boundary conditions to relate the values of 
G at both sides of an interface. For a long time, the 
only boundary conditions for isotropic Green's functions 
(BCIGF) available were spin-independent BCIGF de- 
rived for a S/N tunnel contact—. The only interfacial pa- 
rameter involved in these BCIGF was the tunnel conduc- 
tance Gt of the contact. Such a description is incomplete 
for a general diffusive spin-dependent interface. Spin- 
dependent S/F boundary conditions have been first in- 
troduced in the ballistic regime^—. Recently, many ref- 
erences have used spin-dependent BCIGF— ~- 2 to study 
the behavior of hybrid circuits enclosing BCS supercon- 
ductors, ferromagnetic insulators, ferromagnets, and nor- 
mal metals. These BCIGF, which have been first intro- 
duced in Ref. [H include the G T term of Ref. [IB They 
furthermore take into account the spin-polarization of the 
interface tunnel probabilities through a Gmr term, and 
the spin-dependence of interfacial scattering phase shifts 
through G<p terms. It has been shown that the Gmr 
and Gij, terms lead to a rich variety of effects. First, 
the Gff, terms can produce effective Zeeman fields inside 
thin superconducting or normal metal layers 2 ^ - — , an ef- 
fect which could be used e.g. to implement an absolute 
spin-valve^. In thick superconducting layers, this ef- 
fect is replaced by spin-dependent resonances occurring 
at the edges of the layers 2 ^ Secondly, the G^ terms can 
shift the spatial oscillations of the superconducting order 
parameter in ferromagnets 2 ^— . Thirdly, the G^ term 
can produce superconducting correlations between equal 
spins, e.g. in a circuit enclosing a BCS superconductor 
and several ferromagnetic insulators magnetized in non- 
collinear directions 3 ^. The Gmr terms have been taken 
into account for a chaotic cavity connected to a super- 
conductor and several ferromagnet s 28 ' 29 . In this system, 
crossed Andreev reflections and direct electron transfers 
are responsible for nonlocal transport properties. The 
ratio between these two kinds of processes, which deter- 
mines e.g. the sign of the nonlocal conductanc e 33 ' 34 , can 
be controlled through the relative orientation of the fer- 
romagnets magnetizations. 

In this article, we present a detailed derivation of the 
spin-dependent BCIGF based on a scattering descrip- 
tion of interfaces. Our results thus provide a micro- 
scopic basis for all future investigations of ferromagnet- 
superconductor diffusive heterostructures taking into ac- 
count the spin-dependent interface scattering. To make 
the BCIGF comprehensive and of practical value, we 
make a series of sequential assumptions, starting from 
very general to more and more restrictive hypotheses. In 
a first part, we assume that the contact is fully metallic, 



i.e. it connects two conductors which can be supercon- 
ductors, ferromagnets or normal metals. We consider fer- 
romagnets with exchange fields much smaller than their 
Fermi energies, as required for the applicability of the 
quasiclassical isotropic description. We assume that the 
contact nevertheless produces a spin-dependent scatter- 
ing due to a spin-dependent interfacial barrier Vb- In 
this case, we establish general BCIGF which require the 
knowledge of the full contact scattering matrix. Then, we 
assume that the contact locally conserves the transverse 
channel index (specular hypothesis) and spins collinear to 
the contact magnetization. In the tunnel limit, assuming 
Vb is weakly spin-dependent, we find that the BCIGF in- 
volve the Gt, Gmr, and G^ terms used in Refs. H3-H2, 
plus additional G x terms which are usually disregarded. 
In a second part, we study a specular contact connecting 
a metal to a ferromagnetic insulator (FI) . If we assume 
a weakly spin-dependent interface scattering, we obtain 
the BCIGF used in Refs. [HU We also present BCIGF 
valid beyond this approximation. Note that the various 
BCIGF presented in this article can be applied to non- 
collinear geometries. 

Most of the literature on superconducting hybrid cir- 
cuits uses a spatially continuous description, i.e., in each 
conductor, the spatial dependence of the Green's function 
G is explicitly taken into account. The BCIGF presented 
in this article can also be used in the alternative approach 
of the so-called circuit theory. This approach is a system- 
atic method to describe multiterminal hybrid structures, 
in order to calculate average transport properties 3 ^— 
but also current statistic s 38 ' 39 . It relies on the mapping 
of a real geometry onto a topologically equivalent circuit 
represented by finite elements. The circuit is split up into 
reservoirs (voltage sources), connectors (contacts, inter- 
faces) and nodes (small electrodes) in analogy to classical 
electric circuits. Each reservoir or node is characterized 
by an isotropic Green's function G without spatial de- 
pendence, which plays the role of a generalized potential. 
One can define matrix currents, which contain informa- 
tion on the flows of charge, spin, and electron/hole coher- 
ence in the circuit. Circuit theory requires that the sum 
of all matrix currents flowing from the connectors into a 
node is balanced by a "leakage" current which accounts 
for the non-conservation of electron/hole coherence and 
spin currents in the node. This can be seen as a gen- 
eralized Kirchhoff's rule and completely determines all 
the properties of the circuit. So far, circuit theory has 
been developed separately for F/N— and S/N circuits 3 ^. 
Throughout this article, we express the BCIGF in terms 
of matrix currents. Our work thus allows a straightfor- 
ward generalization of circuit theory to the case of multi- 
terminal circuits which enclose superconductors, normal 
metals, ferromagnets and ferromagnetic insulators, in a 
possibly noncollinear geometry. 

This article is organized as follows. We first consider 
the case of a metallic contact, i.e. a contact between two 
conductors. Section [TT| defines the general and isotropic 
Green's functions G and G used in the standard descrip- 
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tion of hybrid circuits encompassing BCS superconduc- 
tors. Section Hill introduces the ballistic Green's function 
g, which we use in our derivation. Section llVI discusses 
the scattering description of the contact with a trans- 
fer matrix M. Although we consider the diffusive limit, 
the scattering description is relevant for distances to the 
contact shorter than the elastic mean free path. On this 
scale, one can use M to relate the left and right bal- 
listic Green's functions §l and gn. Section [V] presents 
an isotropization scheme which accounts for impurity 
scattering and leads to the isotropic Green's functions 
Gl(r) away from the contact. Section PVTI establishes the 
general metallic BCIGF which relate Gl, Gr and M. 
Section I VIII gives more transparent expressions of these 
BCIGF in various limits. Section [Villi addresses the case 
of a contact with a FI side, in analogy with the treat- 
ment realized in the metallic case. Section llXl concludes . 
Appendix [A] discusses the structure of the transfer ma- 
trix M and Appendix [B] gives details on the calculation 
of the matrix current. Appendix [C] relates our BCIGF 
to the equations previously obtained in the normal-state 
limi t 35 ' 36 . Appendix [D] discusses the BCIGF obeyed by 
the retarded parts of Gum in the collinear case. For 
completeness, Appendix [E| presents the Usadel equations 
in our conventions. 



II. 



GENERAL AND ISOTROPIC GREEN'S 
FUNCTIONS 



From section [TT| to IVIII we consider a planar metallic 
contact between two diffusive conductors noted L (left 
conductor) and R (right conductor) (see Fig. [TTJ) . The 
conductor L[R] can exhibit spin and/or superconduct- 
ing correlations, due to its superconducting order pa- 
rameter A or exchange field E ex , or due to the prox- 
imity effect with other conductors. For the primary 
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description of electronic correlations in L and i?, one 
can use a general Green's function G defined in the 
Keldysh(g)Nambu(8)Spin(g)Coordinate space. In the sta- 
tionary case, G can be defined as 



Gir.F'.:)= I jG(f,f',t 



t ') exp { t£ L_L} (i) 



withiP. 

G{r,r',t-t') = 



G r (r,rV-i') G K {r,r',t-t') 
G a (r,f',t-t') 



(2) 



G r ^(r,r',t-t') = Ti0(±(t~ti))f 3 ({*(t,r),&(ti,r>)}) 

(3) 

and 

G K (r, r>, t-t') = -if 3 ( [*(t, f), &(t, ?>)] ) ( 4 ) 

Here, [.., ..] and {.., ..} denote commutators and anticom- 
mutators respectively, r, f space coordinates, t, t' time 
coordinates, and e the energy. We use a spinor represen- 
tation of the fcrmion operators, i.e. 

&(t,r) = (*}(*, 0, -*!(«, f),*t(«, (5) 

in the Nambu<g)Spin space. We denote by 73 the third 
Nambu Pauli matrix, i.e. T3 = diag(l, 1, —I, —I) in the 
Nambu<g>Spin space. For later use, we also define the 
third spin Pauli matrix i.e. &z — diag(l, — 1, 1, — 1). 
With the above conventions, the Green's function G fol- 
lows the Gorkov equations: 



{ef 3 -H(f)+iA(z)-I 
and 



(z))G(r,?',e) = 5(r,f') (6) 



-'■imp 



(z')) = 6(r,f') (7) 



Here, A corresponds to the gap matrix associated to a 
BCS superconductor (see definition in Appendix |E|) . The 
Hamiltonian H(r) can be decomposed as 



H(r) = H l (z)+H t {p) + V b (z,p) 



(8) 



with z and p the longitudinal and transverse components 
of r. The part H^z) = -{h 2 /2m)d 2 / d 2 z - E ex {z)a z - 
Ep(z) includes a ferromagnetic exchange field E ex (z) in 
the direction Z, and the Fermi energy Ep(z 0) — 
Ef,l(r), whereas the part H t (p) — —{h 2 /2m)d 2 /d'~ + 
V c (p) includes a lateral confinement potential V c (p). 
The potential barrier Vb(z,p) describes a possibly spin- 
dependent and non-specular interface. It is finite in the 
area z G [— &l, bu] only. In the Born approximation, the 
impurity self-energy at side Q € {L, R} of the interface 
can be expressed as T,i mp (z, e) = —ihG(z,e)/2TQ. Here, 
the impurity elastic scattering time tq in material Q can 
be considered as spin-independent due to E ex -C Ep. 
The Green's function G(z,e), which has already been 
mentioned in section |T1 corresponds to the quasiclassical 
and isotropic average of G inside conductor L(R). It can 
be calculated as^i 



G{z,e) = iG(r = R,f' = R,e)/irv Q 



(9) 



with z the longitudinal component of R and vq the den- 
sity of states per spin direction and unit volume for free 
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electrons. Note that we consider geometries where G, 
Tiimp and A are independent of p. 

In this article, we consider the diffusive (i.e. quasiclas- 
sical and isotropic) limit, i.e. 



Structure allowed in the subspaces of: 



E ex ,\A\,e,k B T H/tq < E F 



(10) 



where T is the temperature and k B the Boltzmann con- 
stant. In this regime, the spatial evolution of G(z,e) in- 
side L and R is described by the Usadel equations which 
follow from Eqs. @ and © [see Appendix |E]. The 
characteristic distances occurring in the Usadel equa- 
tions are £ F = (hD F /£ ex ) 1/2 , &v = (hD N /2k B T) 1 / 2 
and £s = (hDg/2 A) 1 / 2 for a ferromagnet F, a nor- 
mal metal N and a superconductor S, respectively, with 
Dq = v F qTq/3 and v f> q the diffusion constant and 



Fermi velocities of material Q. According to Eq. (JTt 
the scale £q is much larger than the elastic mean free 
path = v Fy QTQ. Importantly, the Usadel equations 
alone are not sufficient to describe the behavior of diffu- 
sive hybrid circuits. One also needs to relate the values of 
G at both sides of an L/R interface with BCIGF, which 
we derive in the next sections. 

For the sake of concreteness, we give typical order of 
magnitudes for the different lenghtscales involved in the 
problem. These lenghtscales strongly depend on the de- 
tailed composition and structure of the materials and in- 
terfaces considered, so that the applicability of the qua- 
siclassical isotropic description has to be checked in each 
case. The value of &£(.r) can strongly vary from a few 
atomic layers to a few nanometers if the two materials 
constituting the interface interdiffuse 42 . The mean free 
path, which strongly depends on the impurity concen- 
tration, can be of the order of a few nanometers^. The 
superconducting lenghtscale £s is usually of the order of 
10 nm for Niobiu m 43 i 44 . The Cooper pair penetration 
length can reach £p ~ 10 nm for a diluted magnetic al- 
low like CuNi&, or £jv ~ 1000 nm for a normal metal like 
Cu at T = 20 mK^. 

It is worth to note, at this point, that the derivation 
presented below is not restricted to stationary problems 
on superconducting heterostructures. Actually most of 
the derivations made below do not rely on the specific 
Keldysh structure introduced in Eqs. (HJ-flU) and our re- 
sults can be directly used to describe full counting statis- 
tics in the extended Keldysh technique^ or multiple An- 
dreev reflections 46 . In fact, boundary conditions for arbi- 
trary time-dependent scattering problems have been re- 
cently formulated in a similar spirit^.. However, having 
in mind the many concrete applications of the bound- 
ary conditions in superconducting heterostructures and 
keeping the notation as simple as possible, we derive 
the BCIGF below in the framework of the stationary 
Keldysh-Nambu Greens functions. 
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TABLE I: Meaning of the various decorations used in this 
article for operators defined in the £ space. Each decoration 
corresponds to a particular structure in the transverse chan- 
nels (index n), propagation direction (index s), spin (index 
a), Nambu (index v) and Keldysh subspaces. 



III. BALLISTIC GREEN'S FUNCTION 

Considering the structure of Eqs. ©-((H), for z,z' < 
-hi, or z, z' > bji, one can expand G in transverse modes 



G 1 ' 



(r,r',e) 



y„s,n>s> ( 2 > 2 > £ ); 



Xn(p)Xm(p') 



2Trhy / v n (z,e)v m (z' ,s) 
x exp [isfc„(z)z — isk m (z)z]) (11) 

In this section, we use spin indices a, a' G {j, J,} which 
correspond to spin directions parallel or antiparallel to 
the direction and Nambu indices v' £ {e, h} for 
electron and hole states. The indices s, s' G {+, — } ac- 
count for the longitudinal direction of propagation (we 
use a, a' = ±1, v, v 1 = ±1, and s, s' — ±1 in math- 
ematical expressions). We introduce the wavefunction 
Xn(p) f° r t ne transverse channel n, i.e. H t (p)Xn(p) = 
E n Xn(p), and the corresponding longitudinal momen- 
tum and velocity, i.e. k n (z) = (2m(E F (z) — E n )) 1 / 2 /% 
and v n (z) = hk n (z)/m. Importantly, we have dis- 
regarded the dependences of k n and v n on E ex and 
e due to Eq. [TU] The ~ decoration denotes that the 
Green's function Q can have a general structure in 
the Keldysh<gNambu®Spin(g>Channel(g>Direction space, 
noted £ in the following. In contrast, "denotes the fact 
that Gi(/{) has no structure in the Channel<g)Direction 
sub-space, noted £_ in the following (see the summary of 
notations in Table|T]). Due to Eqs. © and 0, Q(z, z',e) 
is not continuous at z = z /37 i 48 . One can us o 49 i 50 

§(z, z, e) = -in (g{z, z', e) + S 3 sign(z - z')) (12) 

with E 3 the third Pauli matrix in the direction of prop- 
agation space, i.e. (^z)n^ ms > = s5 ss >8 nm 5aa>5vu>1-K- 
Equation (112[) involves a ballistic Green's junction 
g(z,z',e) £ £ which is continuous at z = z' . We will 
see below that this quantity plays a major role in the 
derivation of the BCIGF. 
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For later use, we now derive the equations of evolution 
followed by Q. Inserting Eq. (fTTj) into Eqs. (HO©, one 
can check that, for z ^ z' and z, z' < —b^ (Q = L) or 
z,z' > b R (Q = R), Q follows the equations 



. _ d - 
ihT, 3 VQ— +iA- T, imp (z) 
oz 



g(z,z',e)=0 (13) 



and 

Q(z,z',e) 



- - d - 
-ihS 3 v Q — + i.A - T, imp (z) 



= 0. (14) 



We have introduced above a velocity matrix vq 
with a structure in the channels subspace only, i.e. 

(P L{R)Yns,n' s' = V n( Z ^ 0)S SS > S nn ' S^' 8 VU > Ik, with 1 K 

the identity matrix in the Keldysh space. We have fur- 
thermore assumed that the so-called envelope function Q 
varies smoothly on the scale of the Fermi wave length, 
in order to neglect terms proportional to d 2 /dz 2 and 
d 2 /dz' 2 in Eqs. (H3J[H|)i£. 



IV. SCATTERING DESCRIPTION OF A 
METALLIC CONTACT 



We now define, at both sides of the barrier Vb, two 
ballistic zones (with no impurity scattering) located at 
z e [-c L ,-b L ] and z £ [c R , b R ], with c L(R) - b 



(grey areas in Fig. UU) - In the region z € [— cl,c r ], 
we can disregard the superconducting gap matrix A since 
cl +cr <C £s- Therefore, the electron and hole dynamics 
can be described with the Schrodinger equation 



L(R) 



< 



[sf 3 -H{r)]cf>{?,e)=Q 

or, equivalcntly, 

tf{r,e) [era - H (f)] = 
whose solution has the form 51 



6- CT (r,e) = ]TC : :(z,e) 



Xn(p) 



y/2Trhv n (z) 



(15) 
(16) 

(17) 



in the ballistic zones. Here, tp(r, e) is a vector in the 
Spin<X>Nambu<8>Kcldysh space, and ip(z,e) is a vector in 
the £ space. The index s corresponds again to the lon- 
gitudinal direction of propagation. We have introduced 
v indices in the exponential factors of Eq. pT|) because, 
for the same sign of wavevector, electrons and holes go 
in opposite directions. Therefore, in Eq. (fl7|) . s = +/ — 1 
systematically denotes the right/left going states. One 
can introduce a transfer matrix M.{e) £ £ such that 
-0(c_R,e) = M{e)^p{— cl, e). The matrix M and the 
Landauer-Buttiker scattering matrix can be considered 
as equivalent descriptions of a contact, provided one in- 
troduces small but finite transmission coefficients to reg- 
ularize SA in case of perfectly reflecting channels. This 



regularization procedure does not affect practical calcu- 
lations as illustrated in Section VII. D. Since H{r) does 
not couple electron and holes, M. has the structure 



M = 



M e {e) 
M e (-e)* 



(18) 



in the Nambu subspace. Moreover, At is proportional to 
the identity in the Keldysh space, like H(r). For later 
use, we point out that flux conservation leads to 5 ^ 



7W t E 3 M = M ± 3 M 



(19) 



We now connect the above scattering approach with 
the Green's function description^. With the assump- 
tions done in this section, Eqs. (j6j[7]) give, for z € 
[-c L ,c R ] and z' ^ z 



and 



(ef 3 -H(r))G(r,r',s) = 



(r,r',s)(et 3 -H(r'))=0 



(20) 



(21) 



We recall that in the ballistic zones, G takes the form 
(|llj) . In the domain z' > z, a comparison between 
Eqs. (|20H2ip and (|15H16|) gives, in terms of the decom- 
positions (fTTj) and (JTTJ) 

UG(c a ,c a + 0-,e)U 
= M(s)Ug(-c L +0 + ,-c L ,e)UM\e) 

We have introduced above the transformation U = (1 + 
f 3 +[l— T3]£i)/2 to compensate the fact that the v indices 
do not occur in the exponential terms of Eq. (fTT]) . Using 
Eq. ([T2j). we obtain 



~g R = Mg L Mi (22) 
with g L = g(-c L ,-c L ,e), g R = g(c R ,c R ,e), and 



M 



M e {e) 

EiX e (-e)*Si 



(23) 



in the Nambu subspace. Note that due to Eq. (|19[) . one 

has 



M f t 3 M = M to M f = S a 



(24) 



We now discuss how spin-dependences arise in our 
problem. Due to the hypotheses required to reach the 
diffusive limit [see Eq. (flQ|) ]. we have neglected the de- 
pendence of k n and v n on the exchange field E ex and the 
energy e. Accordingly, we have to disregard the depen- 
dence of M on E ex and e. This does not forbid that M 
depends on spin. Indeed, in the general case, when an 
interface involves a F material which is ferromagnetic in 
the bulk, the transfer matrix M can depend on spin for 
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two reasons: first, the wavectors of the electrons scat- 
tered by the barrier can depend on spin due to E ex , and 
second, the interface barrier potential Vb can itself be 
spin-dependent. Importantly, one can check that E ex 
and Vb occur independently in Eqs. ([6][7]). The value 
of E ex and the spin-dependence of Vb are not directly 
related, because the second depend on properties like in- 
terfacial disorder or discontinuities in the electronic band 
structure, which do not influence E ex far from the inter- 
face. Therefore, nothing forbids to have simultaneously 
E ex <C E F (this can occur e.g. in a diluted ferromag- 
netic alloy like PdNi) and a spin-dependent M, due to 
a spin-dependent interface potential Vb- It is even pos- 
sible to obtain this situation artificially, by fabricating 
e.g. a contact with a very thin FI barrier separating two 
normal metals or superconductors. Note that in spite 
of E ex <C Ep, the exchange field E ex can play a ma- 
jor role in diffusive hybrid circuits by modifying drasti- 
cally the spatial evolution of the isotropic Green's func- 
tion Gp(z,e) inside a ferromagnetic metal F on the scale 
£f [see Appendix [E]. 



V. ISOTROPIZATION SCHEME 

In this section, we show that the Green's function 
g(z, z' — z, e) becomes isotropic in momentum space (i.e. 
proportional to the identity in the £_ subspace) due to 
impurity scattering, when moving further away from the 
contact. One can consider that this process occurs in 
"isotropization zones" with a size dz,(.R) of the order of 

a few £e for side L{R) of the contact (dotted ar- 
eas in Fig. HI]) . Beyond the isotropization zones, quasi- 
particles reach diffusive zones (purple areas in Fig. HI)) 
characterized by isotropic Green's functions G{z, e) with 
no structure in the £_ subspace. We show below that 
g(z,z' — z,e) tends to G(z — z fdL(R),£) at the external 
borders z — ^d,L(R) of the isotropization zones. Note 
that the results presented in this section do not depend 
on the details of the isotropization mechanism. 

We study the spatial evolution of Q in the isotropiza- 
tion zones located at z e [— di,, — cj,] and z £ [cr, cIr], us- 
ing Eqs. (ITBl and (fl4")) . The superconducting gap matrix 
A can be neglected from these Eqs. due to duR) -C £s- 
We thus obtain, for the isotropization zone of side Q and 
z^z' 



d_ 

dz 



2tq 



and 



g{z,z',e) ® ( - s 3Wq^7 



d G(z',e) 



2tq 



= 



(26) 



Due to £q 3> if , one can disregard the space-dependence 
of G(z, e) in the above equations. We will thus replace 
G(z,e) by its value Gq at the beginning of the diffusive 



zone Q, i.e. Gl(A) = G(z = TdL(R),s)- For later use, we 
recall that Gl and Gr. fulfill the normalization condition 



G\ — G R — 1 



(27) 



with 1 the identity in the E space. In the isotropization 
zone of side Q, Eqs. JTJ) and (|25H27p give 



g{z,z',e) = -inP Q [\ Q (z)} 
x 9Q + sign(z - 
with X L ( R )(z) =z± c L{R) and 



(28) 



P q [z] = ch [z/2vqt q ] - S 3 G Q sh [z/2v q t q ] (29) 

for Q 6 {L, R}. Note that the choice of the coordinate 
dh(R) i n Fig-1 is somewhat arbitrary, i.e. defined only up 

to an uncertainty of the order of ie , because there is a 
smooth transition between the isotropization and diffu- 
sive zones of the contact. As a result, Q must tend contin- 
uously to its limit value Qdiff in the diffusive zones. The 
function Gdi ff(z, z' , e) must vanish for \z — z'\ ^> if (see 
e.g. Ref. |54I ). This imposes to cancel the "exponentially 
divergent" terms in Eq. (1281) , which requires^ 



(t 8 + G L ) (g L - S 3 ) = 



(h + E 3 ) ( 



S 3 - G L = , 



(S 3 - Gr) (g R + £ 3 ) = 

(g R - E 3 ) (s 3 + Gr) = 

For z —¥ z fdL(R) we obtain from Eqs. 
finally approaches 



(30) 



(31) 



(32) 



(33) 



that Q 



•g(z,z',s) = (25) | 



Gdiff(z,z',e) = -i7rexp(--=r^ — ) 

2vl(r)1~l(r) 

x (G i(fl) +sign(z-z')S 3 ) , (34) 

so that 9l(r)(z,z' = z,e) tends to G^rm. As required, 

the expression ([3~4")> of Qdiff does not depend on the 
exact choice of the coordinate d^^ and vanishes for 
z - z'\ > if. Equations (|28I33|1 indicate that the de- 
cay length for the isotropization of Ql{r){ z t z ' = z t £ ) 

is max„[(2r7i[£; F:L (7 ? ,) - S„]) 1 / 2 T L(fl) ] = ic^ R \ as antici- 
pated above. Moreover, inserting Eq. (|3~4")l into Eq^. (fTT)) 
leads to an expression of G whose semiclassical and 
isotropic average corresponds to Gur\, as expected^. 
Importantly, from Eqs. (|28H33|) , one sees explicitly that 
Q is smooth on a scale of the Fermi wave length, which 
justifies a posteriori the use of the approximated Eqs. 
(fT3")) and (fT4")) in this section. 
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VI. MATRIX CURRENT AND GENERAL 
BOUNDARY CONDITIONS 



current up to the beginning z — T^Lfft) of the diffusive 
zones allows to identify these expressions with 



Our purpose is to establish a relation between Gl and 
Gr. To complete this task, it is convenient to introduce 
the matrix current^! 



I(z,e) = 

7TTO 



(/ | d_ _ d_ 

' dz dz' 



' . (35) 

This quantity characterizes the transport properties of 
the circuit for coordinate z and energy e. It contains in- 
formation on the charge current (see section [VII Eft but 
also on the flows of spins and electron-hole coherence. 
Note that in this article, e denotes the absolute value of 
the electron charge. Using Eq. (fl"T) and the orthonormal- 
ization of the transverse wave functions Xni the matrix 
current is written as 



I(z,s) = 2iG g Tr„ >s E 3 £(z,z,e) /tt 



(36) 



for z < — or z > b^. Here Tr^^ denotes the trace 
in the £_ sub-space and G q = e 2 /2irh is the conductance 
quantum. Inside the isotropization zones, using Eq. (|28[). 
one obtains 5 - - 

I(z,e) = 2G q Ti niS [t 3 P Q [X Q (z)}~g Q P Q [~X Q (z)]\ . (37) 

Considering that Pq (z) has a structure in the £_ sub-space 
only, and that Pq[~\q(z)]Y, s Pq[Xq(z)} = S3, one finds 



I{z,e) = 2G 9 Tr„„ 



h(R){e) (38) 



at any point in the left (right) isotropization zone. We 
conclude that, quite generally, the matrix current is con- 
served inside each isotropization zone. We will see in 
next paragraph that this property is crucial to derive the 
BCIGF. 

In order to express cjl in terms of Gl and Gr and M, 
we multiply Eq. ([3D) by G L from the left and Eq. ([32) 
by G L M t from the left and by (Aft) -1 from the right. 
Then, we add up the two resulting equations after sim- 
plifications based on Eqs. ([22). ([24) . and ([27). This leads 
to 



He) 



2G q 'TY n 



2D1 1 (G L E 3 + l) - 1 



(39) 



with Dl = 1 + GlM^GrM . A similar calculation leads 
to 



I R (e) = 2G q Tv n 



2D- 1 (g r £ 3 - l) + 1 



(40) 



with D R = 1 + G R {M^)- 1 G L M~ 1 . Equations (Jp) and 
PQ) represent the most general expression for I L ^{e) 

in terms of the isotropic Green's functions Gl^r) and 
the transfer matrix M . The conservation of the matrix 



h{R){s) 



A 



Pl[r\ 



G(z,e) 



dG(z,s) 



dz 



(41) 



z=^d L 



Here, Pl(r) denotes the resistivity of conductor L(R) and 
A the junction area. Formally speaking, Eqs. ([39). ([40) 
and (|4"T) complete our task of finding the general BCIGF 
for spin-dependent and diffusive metallic interfaces. We 
recall that to derive these equations, we have assumed 
a weak exchange field in ferromagnets (E ex <C Ep) 1 as 
required to reach the diffusive limit [see Eq. ([TU) ]. How- 
ever, we have made no restriction on the structure of 
the contact transfer matrix M. In particular, M can be 
arbitrarily spin-polarized, and it is not necessarily spin- 
conserving or channel-conserving. However, at this stage, 
a concrete calculation requires the knowledge of the full 
M (or equivalently the full scattering matrix). Usually 
this information is not available for realistic interfaces 
and one has to reduce Eqs. ([39) - ([4H) to simple expres- 
sions, using some simplifying assumptions. For a spin- 
independent tunnel interface, Eqs. (J3S|)-(jin) can be ex- 
pressed in terms of the contact tunnel conductance Gt 
only, which is a formidable simplification 17 . Another 
possibility is to disregard superconducting correlations. 
In this case, Eq. ([39") and (|4"01 lead to the normal-state 
BCIGF introduced in Refs. I35ll36l (see appendix [C] for 
details). The normal-state BCIGF involve the conduc- 
tance Gt but also a coefficient Gmr which accounts for 
the spin-dependence of the contact scattering probabili- 
ties, and the transmission and reflection mixing conduc- 
tances G l mix and G^^' r which account for spin- torque 
effects and interfacial effective fields 56 . We will show be- 
low that for a circuit enclosing superconducting elements, 
the BCIGF can also be simplified in various limits. 

Note that since the transition between the ballistic, 
isotropization and diffusive zones is smooth, the choice 
of the coordinates and in Fig.l is somewhat 

arbitrary, i.e. defined only up to an uncertainty of the or- 
der of te or a fraction of t% respectively. However, 
one can check that this choice does not affect the BCIGF. 
First, a change of and cr by quantities Scl and Scr 
of the order of a fraction of requires to replace the 
matrix M appearing in Eqs. Q39H41P by Ar M Al, where 

the matrices Ar and Al have a non-trivial (i.e. diagonal) 
structure in the £_ subspace only, with diagonal elements 

Al,u,s = exp[i s 5cl k n ) and Ar^^ s = exp[i s Scr k n ]. 

Since G^(jj) commutes with ^4_r[l] , this leaves the BCIGF 
unchanged. Second, due to Eqs. (|30H33I) . the BCIGF do 
not depend either on the exact values of cLl and cLr. 
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VII. CASE OF A WEAKLY SPIN-DEPENDENT 
S/F CONTACT 

A. Perturbation scheme 

In the next sections, we assume that the transverse 
channel index n and the spin index a =t, i corresponding 
to spin components along Z are conserved when electrons 
are scattered by the potential barrier Vb between the two 
ballistic zones (we use for instance Vb(z,p) = Vo(z)&o + 
Vi(z)az)- In this case, one can describe the scattering 
properties of the barrier with parameters T n , P n , <fn^ R \ 
and dtfn^ defined from 



and 



l L(R), 



= T n (l + aP n ) 



(42) 



and 



rg(r L{R ),n*) = + <?(d^ R} /2) (43) 



with ii(fl) incr the transmission amplitude from side L(R) 
to side R(L) of the barrier and r L ^ na the reflection 
amplitude at side L(R). The parameter P n corresponds 
to the spin-polarization of the transmission probability 

|tfl(i),no-| • The parameters dtp!^ an d dip R characterize 
the Spin Dependence of Interfacial Phase Shifts (SDIPS), 
also called in other references spin mixing angled—. In 
our model, P n and dip n can be finite due to the spin- 
dependent interface potential Due to flux conserva- 
tion and spin conservation along Z, the parameters T n , 

Pn, fn ) an d dipn are sufficient to determine the 
value of the whole Ai e matrix (see Appendix [5] for de- 
tails). Then, using Eq. (|2"3"|) . one can obtain an expres- 
sion for M. We will work below at first order in P n and 
dtfn^. In this case, M can be decomposed as 



M = M°(1 + 6X) 



(44) 



The n th diagonal element of M° in the transverse chan- 
nel subspace has the form, in the propagation direction 
subspace, 



O"0 



(45) 



with R n = 1 — T n . Accordingly, the matrix SX is, in the 
propagation direction subspace, 



SX = 



sx ++ 

SX* 



5X+- 
-8X+. 



with 



SX. 



n+,n+ 



1<TZ 

47^ 



[T n dtfl + (2 - T n )dtp R ) 



(46) 



(47) 



SX n+tn - 



Pn 
\ZRr-i 



\/Rn 



dl Pr. 



(48) 



One can check that Eqs. (|44l) - ([4"5)) are consistent with 
Eq. (|24|). Due to Eq. {23]), the matrices M° and SX 
are proportional to the identity in the Nambu subspace. 
The matrix M° is determined by the parameters T n and 
tfin^ . It has a structure in the £ subspace only. In 
contrast, SX is a first order term in P n and dipn^ R \ with 
a structure in the £_ sub-space but also in the spin sub- 
space. We conclude that the matrices M° and G^m 
commute with each other, whereas SX commutes neither 
with M° nor with Gl(r)- 

We want to express the matrix current of the 
isotropization zones as 



(i) 



L{R) 



(49) 



with (s) and /[Vm (e) zeroth and first order terms in 
SX, respectively. We will mainly focus on the calculation 
of /i(e) because the calculation of /r(e) is similar. To 
develop the expression (|3"9")l of /z,(e), one can use: 



D~ r 



J - JSVJ + o(SV 2 ) 



with 



SV = G L G R Q SX + G L SX^Q G R 



Qo = (M°)tM° 



and 



J = (1 + QoGlGr)- 1 
For later use, we note that 

GrGl + Q 1 



J = 



{G l ,Gr} + Q + Q 



(50) 

(51) 
(52) 

(53) 
(54) 



In the next sections, we will substitute Eq. ([50)1 into 

Eq. (l39|) , to express /| (e) and (e) in terms of the 
scattering parameters of the contact. 

B. Zeroth order component of the matrix current 

We first discuss the conservation of the zeroth order 
matrix current across the contact. From Eqs. (|22l) and 



|, one finds J°(e) = 2iG q Tr n , s \E 3 M°g L (Af°)tj /tt. 

Since M° has a structure in the £_ subspace only, the 
cyclic property of the trace Tr„ :S yields ™(e) = ^(e) = 
/°(e). Hence, the matrix current is conserved across the 
contact in the spin-degenerate case. 



9 



We now calculate I^(e). Since M° commutes with 
G L (R), Eq. m gives 

/<°>(e) = 2G g Tr„, s [j (2E 3 G L + 1 - QoG L G^j\ . (55) 



From Eq. (|451) . one finds 



= -2T n " 



(2T - 



1-T 

1)E 



cos(^)Si + sin(^)E s 



and 



Qo 1 — E3Q0E3 



(56) 



(57) 



In Eq. (1561) . the matrices Eq, Ei and E 2 refer to the 
identity, the first and second Pauli matrices in the 
propagation direction subspace, respectively. We use 

i T o)nV,ms' = T n5 ss '5 nm 8 a( y'5 vv >lK- We find that Q + 
Qq has a diagonal structure in the propagation direc- 
tion space. Therefore, using expression (|54|) for J, and 
performing the trace over the channel and propagation 
direction indices, we obtain 



J< 0) (e) = 4G,£ 



T n [Gr, Gl] 



4 + T n ({G L ,G R }-2) 



(58) 



Equation (|58|) corresponds to the expression obtained in 
Ref. for a spin- independent contact^. This expression 
does not involve any scattering phase shift. 



C. First order component of the matrix current 

We now concentrate on the contribution l£) R \(s) to 

the matrix current to first order in 5X. Equations (|39j) 
and ([50) lead to 



l£\e) = -4G q Tr ntS 



J SV J I + E3 G L 



(59) 



with SV given by Eq. (|5"T]>. Using Eqs. (|46M8[) and ([56l 
I57[) , and performing the trace over the transverse channel 
and propagation direction indices (see Appendix [B] for 
details), one finds 

I^(e) = 2G q (4 + T n ({G L , G R } - 2))" 1 



D. Expression of the matrix current in the tunnel 

limit 



We now assume that the contact is a tunnel barrier 
(T„ <C 1), which seems reasonable considering the band 
structure mismatch between most S and F materials. At 
first order in T„, the matrix currents Il(r){^) take the 
very transparent form 

2^i(e) = Gt [Gr,Gl] +Gmr [{<^z,Gr} , Gl] 

+ iG^ [<tz,Gl] +iG^[GrGl^z + ^zGlGr,Gl] 



iG^ [Gr [&z, Gr] , Gl] 



(61) 



and 



2Ir(e) = Gt [Gr,Gl]+Gmr [Gr, {<tz,Gl}] 

— iG R [(Tz, Gr] — iG^[GlGr&z + OzGrGl,Gr] 



iG^ [Gl [&z,Gl] ,Gr] 



(62) 



We have introduced above the conductance parameter: 



,?>'■> 



Gt/Gq - 2) T n 
Gmr /Gq—y T n P n 
G L ^/G g = -2j: n dv 



L(R) 



The values of the coefficients Gt, Gmr, G 



L(R) 



(63) 
(64) 
(65) 
(66) 

and 



G 



L{R) 



are difficult to predict because they depend on 



the detailed microscopic structure of the interface. These 
parameters can in principle be large compared to G q be- 
cause, although the derivation of Eqs. (f6"T|) and (1521 as- 
sumes that T n , P n and difin^ are small, the definitions 
(|63l) - (|66p involve a summation on a numerous number 
of channels. The parameter Gmr can be finite when 



P n 7^ and the parameters G^ and G£ w can be fi- 
nite due to the SDIPS. From Eqs. (H3J) - <HU) , G£ and G^ 

are likely to be small compared to Gt and G^ R \ This is 
why these coefficients were disregarded so far for studying 
the effects of the SDIPS on the superconducting proxim- 
ity effect. In contrast, it is possible to have G^ > Gt 

M(R) 

~<t> 

potential Vy^i. We also note that the hypothesis P„ C 1 
imposes Gmr "C Gt- We have checked that in the 
normal-state limit, Eqs. (|61I62I) agree with the boundary 
conditions introduced in Refs. [35| and [HI provided the 
reflection and transmission mixing conductances G^^' r 



MR) 



x (4T„P„ [{<tz,Gr},G l ] - iSRndcp^ [& Z ,G L ] 

+ iT n (T n dip% + (2 - T n )d(p%) [G R [<jz, Gr] , G L ] as well as G^ 11 ' < G T , using a spin-dependent interface 
-<T„ (T n dcp% + (2 - T n ) dip%) [[a z , G R ] G L ,G L ]) 
x (4 + T n ({G L ,G fi } -2))" 1 (60) 
A comparison between Eqs. ([39]) and ([40]) indicates that 



the expression of l2^(e) can be obtained by multiplying ^^^.^^ ^-mix 

H J,-. T /A\ and G mix appearing m these boundary conditions are re- 



the expression ([60]) of I L 1} (e) by -1, replacing d^ w by placed by their developments at first order in T n , P n , and 



MR) 



difn (L \ and G L(R) by G R[L) 



Note that the expressions 



of 1^ (e) and Zjj (e) involve the SDIPS parameters dip^ 
and difin but not the spin- averaged phases (p^ and tp R . 



d(p n i.e.: 



G L m [^ r ^ (Gt/2) + *(G^/2) + 2zG^> 
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and 



spin-dependent situation, one finds from Eqs. ([22]) and 



G t mix ^(G T /2) + l (G L x +G R ) 



that Ir{e) 



2iG q Tv n 



X 3 Mg L (My /n. Since M 



(see Appendix ICl for details). 

We now briefly review the physical effects of the coeffi- 
cients Gt, Gmr, and . The term in Gt in Eqs. (|6"Tj) 
and corresponds to the term derived in Ref. [l3 for 
superconducting/normal metal interfaces. This term is 
responsible for the superconducting proximity effect oc- 
curring in a normal metal layer or a ferromagnetic layer 
in contact with a superconductor. The parameter Gmr 
accounts for the spin-dependence of the contact tun- 
nel probabilities, and thus leads to magnetoresistance 
effect o 28 ' 29 ' 32 . In a ferromagnet F subject to the proxim- 
ity effect, the ferromagnetic exchange field causes spatial 
oscillations of the isotropic Green's function G, which re- 
sults e.g. in spatial oscillations of the density of states 

L ( R) 

of F. It has been shown that the G, terms can shift 

these oscillations 2 ^—. The G^ R ^ terms also induce Zee- 
man effective fields inside thin superconducting or nor- 
mal metal layers^ 2 - - — . In principle, in non-collinear ge- 
ometries enclosing several ferromagnetic elements with 
non-collinear magnetizations, the SDIPS terms can in- 
duce spin-precession effects. 

Note that, so far, we have considered that the inter- 
face potential Vj, is spin-polarized along the Z direc- 
tion. In the general case, due to interface effects, it 
is possible that the spin-dependent part of the inter- 
face potential Vf, is polarized along a direction m dif- 
ferent from the bulk exchange field direction of con- 
tacts L or R. It is also possible that the contact be- 
longs to a circuit enclosing several ferromagnets with 
non-collinear magnetizations, or ferromagnets with a spa- 
tially dependent magnetization direction. In these cases, 
Eqs. (|39fl41|) are still valid. One can furthermore gen- 
eralize the BCIGF (JHOJ), dHJ) and to an arbitrary 
spin reference frame (<Jx, 0Y , &z) by replacing &z by 
((1 + f 3 )a z (m ■ a)a z + (f 3 - l)ov(m • a)a y ) /2. 

As we have already explained in Sec. IIV1 the use of 
transfer matrices for the derivation of Eqs. (|6H66I) al- 
lows to obtain results for the T n — > limit, which must 
be performed after an explicit calculation of the BCIGF. 
From Eq. (|55|) , even if a channel n is perfectly reflected at 
the L/R boundary, it can contribute to the matrix cur- 
rent due to the spin dependence of the reflection phase 
dipn {R) - We will recover this result in Sec. IVIII Al for a 
S/FI contact, using an approach without transfer ma- 
trices. 



E. Discussion on the matrix current conservation 
and the spin-dependent circuit theory 

In this section, we discuss the non-conservation of the 
matrix current in the general case. We have already seen 
in section [VII Bl that the full matrix current is conserved 
across an interface in the spin degenerate case. In the 



has a structure in the spin subspace, the cyclic property 
of the trace Tr„ iS cannot be used anymore to relate 
and Ir{s). Hence, nothing imposes I hip) = -Tr(e) m the 
general case. Reference 24 illustrates that in the case 
of a simple S/F bilayer with a homogeneous magnetiza- 
tion in F, /i(e) ^ Ir{s) is already possible. Note that 
7^ Ir{ s ) does not violate particle current conserva- 
tion through the interface, although the average current 
flowing at side Q of the contact is determined by Iq, i.e. 

1 f°° 

(Iq) = ^ J deTr va {f 3 l£ (e)} (67) 
Indeed, the above equation leads to 

(II) = ^ <feT WCT {f 3 £ 3 g*} 

and 

C f°° ( ~\ 

(Ir) =-^J deTr nsva {f 3 £ 3 Msf M f ) 

Since M is proportional to the identity in the Keldysh 
space, one can use the cyclic property of the trace Tr nsV(T 
in the above equations, to show that (I if) — (Ir)- It is 
important to point out that the non- conservation of the 
matrix current at the L/R boundary does not affect the 
applicability of Eqs. (|39H4ip . The fact that the matrix 
current is not conserved through a spin-dependent inter- 
face has the obvious reason that only charge conservation 
is required by fundamental laws, whereas other quantities 
are not conserved in general. It depends on the symmetry 
of the Hamiltonian describing the barrier, which quanti- 
ties are conserved in addition to charge. If the barrier 
potential is spin-independent, all elements of the matrix 
current are conserved. In general, this is not the case 
anymore for spin-dependent barriers. An extreme case 
illustrating this situation is provided by an interface be- 
tween a FI and a metal. In the FI, the concept of a 
matrix current does not even exist, although the FI in- 
fluences the adjacent metal due to the proximity effect. 
We will discuss this case in section IVIIII 

The BCIGF derived in this article allow to generalize 
the " circuit theory" of Ref. H3 to the case of multitermi- 
nal circuits which enclose superconductors, normal met- 
als, ferromagnets and ferromagnetic insulators. In the 
approach of circuit theory, a system is split up into reser- 
voirs r (voltage sources), connectors c (contacts, inter- 
faces) and nodes n (small islands) in analogy to classical 
electric circuits. Each reservoir or node is characterized 
with an isotropic Green's function with no space depen- 
dence, which plays the role of a generalized potential. 
Circuit theory requires to apply generalized Kirchhoff's 
rules on the matrix current /. We have seen above that / 
is not conserved through the contacts in the general case, 
but this is not a problem since we know how to express 
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the matrix current at both sides of the contact. We will 
note /" the matrix current flowing from the connector c 
into node n, which is given by Eqs. (|39|) or (|40|) . One must 
be careful to the fact that the matrix current is not con- 
served either inside the nodes due the terms on the right 
hand side of the Usadel Eq. (|E1[) . To compensate for the 
non-conservation of / inside node n, one can introduce a 
leakage matrix current 



jn 

leakage 



4ttG> K [-iefa + A + iE ex a z ,G n ] (68) 



which accounts for the "leakage" of quantities like for in- 
stance electron-hole coherence or spin accumulation. In 
the above expression, G n , A, and E ex , refer to the val- 
ues of the isotropic Green's function, gap matrix, and 
exchange field inside n, and V n is the volume of the 
node. The leakage matrix current I™ eakage can be viewed 
as flowing from an effective " leakage terminal" . It must 
occur in the generalized Kirchhoff 's rule for node n, i.e. 



leakage 



J" = 



with the index c running on all the contacts connected 
to node n. We refer the reader to Refs. [37lf60l for more 
details on circuit theory. 



VIII. CONTACT BETWEEN A METAL AND A 
FERROMAGNETIC INSULATOR 

A. Boundary conditions 

In the case of a contact between a metal and a ferro- 
magnetic insulator, one can perform a calculation similar 
to the one of the metallic case without using the trans- 
fer matrix M but a simpler "pseudo" transfer matrix M 
which involves only reflexion phases against the FI (see 
definition below). This facilitates a developement of the 
BCIGF at higher orders in the SDIPS parameters. We 
assume that the ferromagnetic insulator is located at the 
right side (z > 0) of the contact, and that the left side 
L is a BCS superconductor, a normal metal, or a ferro- 
magnet, which satisfies Eq. (flU]) , We divide L into a 
ballistic zone, an isotropization zone and a diffusive zone 
like in Figure [TTJ We directly consider the case where 
the contact locally conserves the transverse channel in- 
dex and spins along Z. In this case, one can introduce 
a reflection phase shift tp n + adip n /2 such that the left- 
going and right-going quasiparticle wavefunctions in the 



n channel of L are related by 



V#,-(- c £> e ) 



+<y- 



Using this relation, one can check that the calculations 
of sections |TV] to [VI] can be repeated by replacing the 
ballistic Green's function by Eig^Ei, Gf> by Gl, and 
the transfer matrix M by a pseudo transfer matrix M = 
M° exp((5X). The n th diagonal elements of M° and 5X in 



the transverse channel subspace write 

M° = [cos (ip n ) S + i sin (ip n ) ± 3 



and 



SXn = idip n E 3 az/2 ■ 



(69) 



(70) 



Since M° commutes with G L and SX and (M°)tM° = 1, 
we find 



I L (e) = 2G 8 Tr„,.. 
with 

SY= (G L e 



G L E 3 



1-1 



l)/2 



(71) 



(72) 



Hence, quite generally, the spin-averaged reflection 
phases <p n do not contribute to Equation (fTTj) 

can be traced out numerically. Alternatively, one can 
achieve further analytical progress by expanding 
with respect to the spin-dependent part SY. We have 

1 + 5Y) = 1 + E„ {-^Y- Therefore, at fourth 
order in dip n we obtain 

2Ji(e) = iGfai [<tz,Gl] +G^2 [&z,Gl&zGl] 

> 2" 



+ Gfa4 



Gl \pzGh) 
&z, Gl (ctzGl) 



(73) 



with the conductance parameters 

* — ' n * — • n 

£ ^ n £ j n 

G^/G q = Y, d ^n/Z 

t, — J n 
* — 

In the normal-state limit, we have checked that Eq. 
agrees with the BCIGF presented in Refs. and [36 
(see appendix [C] for a detailed comparison) . The term 
in GVi of Eq. (1731 has already been used in Refs. [12 
and |23| . At first order in d<p n , it is the only term con- 
tributing to I i,, and it can be recovered from Eqs. (|61[) 
and ([55 |) -([e>o ]l by using T n = and d^ = d<p n . At 
higher orders in dip n , the value of G^^ is renormalized 
and new terms occur in the expression of Ii . The second 
order term has a straightforward interpretation, since it 
has exactly the same matrix structure as the self-energy 
due to scattering by paramagnetic impurities in a nor- 
mal meta l 61 i 62 , or due to magnetic disorder along the Z 
direction in a ferromagnet 63 . The scattering of Cooper 
pairs at the spin-active interface leads to a coupling be- 
tween spin-singlet and spin-triplet components, which, 
due to the random scattering at second order leads to 
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pair breaking. In a similar fashion, we can understand 
the higher order terms in Eq. [73] as a result of multiple 
scattering at the S/FI interface. Note that in this sec- 
tion, we have assumed that the FI side of the contact is 
magnetized along the Z direction. If the FI is magne- 
tized along a direction m ^ Z, one can describe the con- 
tact in the spin reference frame (ax , <Jy , &z) by replacing 
a z by ((1 -I- f 3 )a z (m ■ a)a z + (?3 - l)a y (m ■ a)a y ) /2 in 
Eq. (ESI). 



B. Example of a S/FI bilayer 

To illustrate some effects of the G^i coefficients, we 
now consider the case of a S/FI bilayer, with S located 
at z £ [0,ds] and FI at z > ds- Throughout this sec- 
tion, we replace the energy — ie appearing in the Usadel 
equation by —ie + T, where the phenomenological colli- 
sion rate T accounts for inelastic processes^. Inside S, 
the retarded part of the isotropic Green's function can 
be parametrized with a so-called pairing angle such 
that 



G r 



cos(A^) 



sin(A^) 





cos(A^) 
sin(Af) 






sin(A^) 
- cos(A^) 




sin(A^) 



— cos( 



,s(Af) 



Let us first assume that ds *C £s> so that one can use 
the quadratic approximation Af (e, x) = A° — /3 CT (x/£s) 2 
and a constant superconducting gap A(x) = Ao inside S 
(see e.g. Ref. For z £ [0, ds], the Usadel equations 
(see appendix IE! lead to 



A cos(A°) + (ie-r)sm(A°) 



2A 



BCS 



(78) 



We have introduced above the bulk BCS gap Abcs of 
S. The value of A 6 can be found by identifying Eq. (J75J) 
with Eq. fl73J), i.e. 

2p a d a / Xs = i70,i°- sin ( A ") + 70,2 sin(2A°) 

+ i^a sin(3A°) + 7<M sin(4A°) (79) 

(see Appendix [D] for details) . We have introduced above 
70,i = G<f>,it;sPs/A. Note that the value of A must be 
calculated self-consistently with A°, see e.g. Ref. l24l We 
will first consider the case G^,2 = G^.3 = G^,4 = 0, for 
which Eqs. ((7511791) yield 



A" = arctan 



A 



-ie + T + nV>,iCTd7 



(80) 
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From the above Eq., G^i induces an effective Zeeman 
field H e ff = 2i704^ s A B cs/d s g^B inside a thin 5 layer, 
like the G^ (fl) terms of section IVIID! ^. The den- 
sity of states (DOS) in the S layer can be calculated G^.4 terms. Assuming A <C 27r, the linearization of Eq. 



FIG. 1: Density of states N(e, x = ds) in a S layer contacted 
to a FI. The black dashed lines correspond to V — 0.1 and 
7c/>,2 = 0, and the blue full lines correspond to T = 0.01 and 
7<#>,2 = 0.03. The left panel corresponds to ds/£s — 0.5 and 
the right panel corresponds to ds/Cs = 3. In all cases, we 
have used 7^1 = 0.15, 7^,3 = 7^,4 = and k B T = 0.1A S cs- 



as N(e,x) = N ^ Re[cos(Af (e, x))]/2, with N the 
normal-state density of states. The G^i-induced effec- 
tive Zeeman field H e f f splits the superconducting peaks 
of the DOS, as shown by the black dashed line in Fig. [TJ 
left panel. Spin-splitting effects in S/F systems were first 
intuited by De Gennes from a generalization of Cooper's 
argumen t 65 ' 66 . Later, Ref. [l^ has confirmed from a qua- 
siclassical approach that the SDIPS can induce a spin- 
splitting of the DOS in a ballistic S/FI bilayer with a 
thin S. However, the effect found by Tokuyasu et al. is 
qualitatively different from ours. Indeed, in the ballistic 
limit, Tokuyasu et al. find that the S/FI bilayer differs 
from a S layer in an external field because the SDIPS 
induced spin-splitting effect depends upon the quasipar- 
ticle trajectory. In contrast, in the diffusive limit, we 
obtain a true effective Zeeman field H e f f which appears 
directly in the spectral functions. On the experimental 
side, spin-splitted DOS were observed in superconduct- 
ing Al layers contacted to different types of FI as soon as 
1986 (see Refs. [67HZ3)- However, the inadequacy of the 
ballistic approach of Tokuyasu et al. for modeling the 
actual experiments was pointed out in Ref. [70l In fact, 
most of the experiments on A\/FI interfaces were inter- 
preted by their authors in terms of a diffusive approach 
with no SDIPS, and an internal Zeeman field added ar- 
bitrarily in the Al layer (see Refs. [69Tj7lh . Our approach 
provides a microscopic justification for the use of such an 
internal field. Remarkably, it was found experimentally 70 
that the internal field appearing in S scales with dj 1 , in 
agreement with our expression of H e ff. 

We now discuss briefly the effects of the G^p, G^ and 
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FIG. 2: Top panel: Density of states N(e,x = ds) in a S 
layer contacted to a FJ. The red dashed line corresponds to 
r = and 70,2 = in both panels. The black full lines in the 
top panel correspond to T / and 7,^2 = and the blue full 
lines in the bottom panel correspond to F = and 7^,2 7^ 0. 
In all cases, we have used ds/Cs — 0.5, 7^,1 = 7^,3 = 7^,4 = 
and k B T = 0.1A flC s. 



([75)) leads to 



An 



£ A i(7*. 1+37*, 3)^+27*, 2+47*, 4 

t,sA-BCS 3- 



— IE 



(81) 

Therefore, in this limit, G<j, 3 contributes to the Zeeman 
effective field like G^,x- Moreover, the coefficients G^,2 
and G0,4 lead to a decoherence effect similar to the de- 
coherence induced by the Y term. However, it is clear 
from Eq. (|79|) that this picture is not valid in the gen- 
eral case. Let us focus on the effect of G^- From (fTTj|). 
in the non-linearized limit, 7^.2 occurs together with a 
sin(2A°) in the expression of f3 a . Therefore, as already 
pointed out in section IVIII A[ in the general case, it is 
more relevant to compare the effect of 7^2 to that of 
paramagnetic impurities which would be diluted inside 
S. The analogy to magnetic disorder can be understood 
as arising due to successive reflections on the S/FI inter- 
face with random spin-dependent phase shifts. To study 
the effect of 7^.2 in the general case, we have calculated 
the density of states N(e, x) numerically. Our code takes 
into account the self-consistency of the superconducting 
gap A(x) in the S layer and is valid for arbitrary val- 
ues of dgH. Figure [2] compares the effect of Y 7^ (top 
panel) with the effect of G^i 7^ (bottom panel), for 
G<f, t i = 0. As expected, we find that the effect of G^ on 
the DOS of a thin 5 is quite similar to the effect of para- 
magnetic impurities which would be diluted inside the 
bulk of ff 73 ' 74 . First, a weak G^.2 widens the BCS peak 
in a way which is qualitatively different from Y since the 



x 



1 

n 

ll 

ll 
ll 
1 1 
1 1 
1 1 

1 1 
1 , 

|A\ 
I V 
1 V 
I \ 


d s /^ s =0.5 

-— d'^0.8 
d s /§ s =1 -5 

— - d s /^ s =8 










/ / /J: x=o 

_ / y _ - 


fl'i 



e / A, 



BCS 



e/A, 



BCS 



FIG. 3: Density of states in the S layer contacted to a FI, 
for x = (left panel) and x = ds (right panel), and different 
values of ds/^s- In all cases, we have used 7,^2 = 0.1, 7^.1 = 
70,3 = 70,4 = 0, T = 0.025 and k B T = O.lAscs. 



curvature of the DOS for e < Abcs has opposite signs 
in the two cases. Second, even a very small Y leads to 
a finite zero-energy DOS, whereas a small G^ j2 reduces 
the gap appearing in the DOS but preserves N(e, x) — 
for small energies. For larger values of G^, we expect a 
gap suppression in the DOS (not shown in Fig. [2j. Note 
that in Fig. [2j for 7 0j2 = 0.1, r = and d s /£s = 0.5, 
the gap of the DOS would be reduced but still finite (not 
shown). In these conditions, using a small Y can trig- 
ger the gap suppression as shown by the black full line 
in Fig. |21 right panel. In the G^i 7^ case, the ef- 
fects of G^,.2 on a thin S remain qualitatively similar, 
in particular, the gap reduction in the DOS still occurs, 
combined with the G^i-induced spin splitting. Figure 
[U left panel compares a case with G§$ = and a large 
r (black dashed line) with a case with G^^ 7^ and a 
small r (blue full line). The two cases can easily be dis- 
criminated due to the different curvatures in the DOS. 
Importantly, the analogy between a paramagnetic impu- 
rity term and the G^ term is not complete since G^ 
occurs in the BCIGF whereas paramagnetic impurities 
would contribute directly to the Usadel equation. This 
discrepancy is revealed by the dependence of the DOS on 
ds- Figure [3] presents the DOS at the left and right side 
of the S layer for different values of ds and G^i = 0. We 
obtain a strong dependence of the G^ DOS-widening 
on x and ds- First, for g?s/£s — 0.5, the DOS at the 
left and the right sides of S (left and right panels) are 
almost identical, with a suppressed gap for the parame- 
ters we consider. When ds increases, the gap reappears 
in the DOS. For d s > <fs, the DOS at the left side of 
S tends to the bulk BCS DOS, with no effect of G^ |2 , 
whereas the DOS at the right side of S still has a reduced 
gap. In this limit, one can check that the reduction of 
the gap occurs for a slab of S of thickness ~ £s near 
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the S/ FI interface. In contrast, paramagnetic impuri- 
ties would affect the bulk of S. Let us now consider the 
case G^i 7^ and ds large. In this case, Ref. |25| has 
shown that the G^i-induced spin splitting of the DOS 
can persist in a slab of S of thickness ~ £s near the S/ FI 
interface. The right panel of Fig. [TJ shows an example 
of DOS at x — ds for ds = 3£s, in the case G^i 7^ 0, 
G(j, t 2 = and a large T (black dashed line), and in the 
case Gc/,,1 7^ 0, G<^2 7^ and a small T (blue full line). 
In the first case, the G^i-induced spin splitting of the 
DOS is not visible anymore because H e jf scales with 
1 /ds and thus becomes too small compared to the large 
value of T used. In the second case, the double gap split- 
ting is still slightly visible as a cusp in the DOS curve 
because the G/,,2 DOS-widening also decreases with ds- 
The effects of the G<^3 and G<^4 terms in the general case 
will be presented elsewhere. Before concluding, we note 
that in circuits enclosing several FI with non-collinear 
magnetizations and BCS superconductors, it has been 
found that the G<^i term can induce spin-precession ef- 
fects which lead to superconducting correlations between 
equal spins^i. 



IX. CONCLUSION 

To model the behavior of electronic hybrid circuits, a 
proper description of the contacts between the different 
materials is crucial. In this article, we have derived gen- 
eral boundary conditions relating isotropic Green's func- 
tions at both sides of the interface between two diffu- 
sive materials [Eqs. flU}, (0DJ and ([ST])]. These BCIGF 
are valid for a circuit enclosing superconductors, normal 
metals, and ferromagnets, in a possibly non-collinear ge- 
ometry. In general, they require the knowledge of the 
full contact scattering matrix, an information usually not 
available for realistic interfaces. However, we have shown 
that in the limit of a specular tunnel contact with weakly 
spin-dependent scattering properties, the BCIGF can be 
expressed in terms of a few parameters, i.e. the tun- 
nel conductance Gt of the contact, a parameter Gmr 
which accounts for the spin-dependence of the contact 

scattering probabilities, and four parameters G^ fl ' and 

L ( R) 

G x which are finite when the contact exhibits a SDIPS 
[Eqs. (|oT]) and (152]) ]. In the case of a contact with a FI 
side, we could simplify the BCIGF for a stronger SDIPS 
[Eq. [73]. We believe that the various spin-dependent 
BCIGF derived in this article represent a solid basis for 
further developments on superconducting hybrid circuits. 
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Appendix A: Scattering description of a specular 
and spin-conserving contact 

1. Structure of the electronic scattering matrix 

In this section, we assume that the transverse channel 
index n and the spin index a e {1\4-}i corresponding to 
spin components along Z, are conserved when electrons 
cross the potential barrier Vf, between the two ballistic 
zones. In this case, the electronic scattering matrix S e is 
diagonal in the (transverse channel)® spin subspace. The 
scattering submatrix associated to electrons with spins a 
of the n th transverse channel writes 



si: 



.no 

tt 



(Al) 



Here, 



r L(R), 



denotes the reflection amplitude 



at side L(R) of the barrier and t 



R(L),na 



the 



transmission amplitude from side R(L) to side 
L(R). Flux conservation imposes, for a e {til} 7 
J2qe{LM} (arg(r Q , n(T ) - &rg(t Q , n(T )) = tt[2it] and 
1 - \TQ, n <r\ = \tQ,na\ 2 = T na . In addition, spin- 
conservation along Z allows to map the scattering 
description of each spin component a onto a spinless 
problem. Time reversal symmetry in each of these 
spinless problems implies arg(t[, )nCT ) = arg(ifl ncr ). 
Therefore, one can use, without any loss of generality 



S^,„ — 



,)/2 



VI - T nr7 e l 



with tpna R ^ = arg(r L ( fl ) nrT ). The matrix S^ a is entirely 
determined by T n<7 , ip^o- an d fna- ^ n this article, we use 

L( R) 

the parametrization T na = T n (l + aP n ) and tp n & = 



L(R) 



a(dtpn iR) /2) [Equations g2][43])]. 



2. Expression of the transfer matrix with 
scattering parameters 

In this section, we assume that the transmission am- 
plitudes tzjm, n <7 are finite. With the hypotheses made 
in section [At] the electronic transfer matrix A4 e is also 
diagonal in the (transverse channel) <g>spin subspace. In 
the propagation direction subspace, the submatrix M. e no 
has elements 



,52 



M 



t 

L,na 



Ml 



(tR,na) 



{tR.na) 1 TL.na 



(A2) 



(A3) 



(A4) 
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Mi 



(t 



R,na) 



(A5) with 



We have used above +/— to denote the right/left-going 
propagation direction. Using Eqs. ([25]) . ([A2HA5P and the 
parametrization introduced in section lA 11 one can obtain 
and expression for the matrix M in terms of T n , P n , tpn 
and d(fin^ ■ At first order in P n and difin^ R \ this leads 
to the expressions (|44H48[) . 



Appendix B: Calculation of I^\e) for a S/F contact 

In this section, we give details on the calculation of the 
contribution I^\s) to the matrix current Il{s) to first 
order in SX. Using Eq. ([55)1. one can rewrite Eq. ([5T)[) 
as 

l£\e) = -4G g Tr n |f + % [{G l ,G r } - 2])"' 
T U + T [{Gl,Gr}-2] 



x Tr, 



W 



(Bl) 



The central term 



W = (GrGl + Qa 1 ) SV (GrG l + Qo x ) (l + S 3 G L 

(B2) 

of this expression can be decomposed as W — 2j=i Wj> 
with 

= QoSXG R G L + GrSX^Q G l 

+ GlGrSXQv 1 + GlQ^SX^Gr (B3) 

W 2 = Qo^XQo 1 + ^^Gfl 

+ G L & R 6XG R 6 L + G L Q^SX^Q Q G L (B4) 



and VU 



3(4) 



W a 1 ( 2 )S3Gl- We now develop the trace over 
the propagation direction index s in Eq. (jBll) . using ex- 
pressions (l46l) . (|56)l and ([57]) . and keeping in mind that 
Gl and (5^ have no structure in the £_ subspace. We find 
Tr s (8X) =Tr s (5X*) ^Ty^QoSXQ- 1 ) =Tr^, W 

Q ) = 0, so that TV S (W 2 ) = 0. Due to Eqs. ([56157]) , we 
find Tr s (Q E 3 ) =Tr s (Q5" 1 E 3 ) = 0. Hence, OT+ + and 
the diagonal elements of Qq and Qo 1 do not contribute 
to Tr s (Wi). In contrast, the development of Tr s (W 3 (4)) 
involves both SX + + and SX^ We finally obtain 



Tr, 



Tr, 



Wx 



TU :i 



[{A,Gr},Gl\ (B5) 



= [[C,G r ]G l ,G l ] (B6) 



Tr., 



W 4 



= 2 [B — Gfl [<5A > ++ , Gfl] , Gl] (B7) 



G = 2Go,++<5^4 



(B8) 



(B9) 



5 = (Qo.++ + Qo,+-Qo,-+) SX++ + Q ,++F - 5X++ 

(BIO) 

Expressing Qq and SX in terms of the scattering param- 
eters T n , P n , (pn (R \ and dipn (R) [see Eqs. gBJUgJ and 
(|56p ] , and developing the trace on transverse channels in 
Eq. ([Bill , we obtain the expression (|60[) for (e). 



Appendix C: General boundary conditions in the 
normal-state limit 



When there are no superconducting correlations in the 
circuit, the isotropic Green's functions Gl(r) write, in 
the Keldysh space: 



G 



L{R) 



T3 K L(R) 
-f 3 



(CI) 



In this limit, the elements D L and D R appearing in the 
general BCIGF (|39I40[) take a simple form. For instance, 



one finds, in the Keldysh space, 



Nl -f 3 N L (M^KrM - K L M^M) N l 
Nl 



with N L = (l + M^M) . A similar expression can be 
obtained for D^ 1 by replacing M by M~ x and K l \r\ 
by Kril] ■ For comparison with sections I VIII and IVIII1 
we specialize to the case of a specular contact con- 
serving spins along the interface magnetization. Equa- 
tions ([391401) give, for the Keldysh electronic component 
of the matrix currents, 



/£ ' e (e) = 2G g Tr„ -t R K R t R + K e L - r L F£r[ 
and 



V ( e ) 



2G g Tr„ 



tLklt{-k R + r R k R r R 



(C2) 



(C3) 



Assuming that the contact is magnetized along Z , we 
obtain 



fK,e 
L L(R) 



{e)/2={[Gt/2]+Gmr)^ 
+ ([G t /2}-Gmr)h 



K 



K 

k 



e,t,T 
R 

R 

e .4-.t>,. 



n 



t GLzx^k^ui t (Glnix) Hk^L^ 



±G 



L(R) 



(C4) 
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with u t (4.) = 1 ± (&z/2), 
and 

We have checked that in the normal state limit, Eq. (|TTj) 
leads to Eq. JC2]) with t R = 0. Equations ([021), (fC3j) 
and (|C4I) are in agreement with the normal-state BCIGF 
presented e.g. in Eq. (2) of Ref. H(| up to a prefactor 
which corresponds to our conventions 75 . Importantly, 
the derivation of these equations requires no particular 
assumptions on the values of t L ^ na and r L ^ niJ . In 
the normal state-limit, a strong spin relaxation is often 
assumed in F, so that the G^^ term is disregarded (see 
e.g. Eq. (5) of Ref. |H). When the circuit includes su- 
perconducting elements, the expressions of DJ^ an d D^ 1 

involve e.g. factors (l + G a L (r) G R {r) M) 1 instead of 

Nl- This is why the superconducting BCIGF are difficult 
to simplify in the general case. 



Appendix D: Equilibrium boundary conditions in 
the case of superconducting correlations between 
opposite spins only 



S V b F 




FIG. 4: Scheme representing two particular types of Andreev 
reflection processes which can occur on a S/F interface mod- 
elled like in figure[Il] The ballistic, isotropization and diffusive 
zones of S and F are represented by grey, dotted and purple 
areas respectively. Full (dashed) lines represent trajectories 
of electrons (holes) from the t (i) spin band. The supercon- 
ducting gap A is taken into account in the diffusive part of S 
only, so that one can consider that Andreev reflections (AR) 
occur at the interface between the diffusive and isotropiza- 
tion zones of S. The upper part of the scheme represents and 
electron incident from the F side, which is transmitted by the 
barrier Vb as an electron, Andreev-reflected on the diffusive 
part of S as a hole, and transmitted again by Vb as a hole, 
before joining the diffusive part of F again. The probability 
associated to this process is proportional to T^(l — P%). The 
lower part of the scheme represents a more complicated tra- 
jectory which also involves two reflections on VJ,. The joint 
probability of these reflections is (1 — T n ) 2 — T^P 2 . 



This appendix presents the boundary conditions 
obeyed by the retarded part of the isotropic Green's func- 
tions, in a case where there are superconducting correla- 
tions between opposite spins only. This situation occurs 
e.g. when all the ferromagnetic elements of the circuit 
are magnetized in collinear directions. For simplicity, we 
assume that no phase gradient is present in the system. 
The conventions chosen in section [TT] give, inside conduc- 
tor Q, 



s(A?) 











cos(A^) sin(A^) 
sin(A^) -cos(Aj) 
sin(A^) 



sin(A^) 



-cos(A^) 



(Dl) 

with A^ = A® in the spin-degenerate case. For a metallic 
contact, using Eqs. (|4"Tj) . (|61l) and (1621) . one obtains 



A dA% 
Pl dz 



G T sin(A^ — A^) + iG^a sin(A^) (D2) 
+ 2i sin(A«)a (G£ - G\ cos(A^ - A*)) 



and 



Interestingly, the Gmr term vanishes from Eqs. (|D2ID3[) . 
so that the tunnel rate polarization P n does not con- 
tribute to the equilibrium BCIGF (we have checked that 
this property remains true when phase gradients occur 
in the system). This result may seem surprising since 
Andreev reflections, which modify the equilibrium den- 
sity of states in a superconducting hybrid system, are 
suppressed when P n is strong 76 . However, one should 
keep in mind that an Andreev reflection process on the 
L/R interface involves together the transmission [or re- 
flection] of an electron and a hole from opposite spin 
bands through the V& barrier (see Figure [3}. These two 
processes have a joint probability which involves P%. In 
contrast, single quasiparticle processes, whose probabil- 
ities involve P n at first order, do not matter at equilib- 
rium. We conclude that P n vanishes from the equilibrium 
BCIGF (|D2ID3[) because we have derived these Eqs. at 
first order in P n . Note that even in this limit, P n docs 
not vanish from the boundary conditions obeyed by the 
Keldysh part of the isotropic Green's functions (see e.g. 
Eq.EU. 

For completeness, we mention that in the case of a 
contact with a FI side, using Eqs. (|4"Tj) and ([73]) . one 
obtains 



A dA* 



G T sm{A L a - A«) - iG$a sin(A«) (D3) 



- 2iasin(A£) (G£ - G£cos(A* - A£)) 



A_dA^_ 

Pl dz 



■■ iG^a sin(A%) + G£ 2 sin(2A£) 
iG^o sin(3A£) + Gl A sin(4A^) (D4) 
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Appendix E: Usadel Equations 

For completeness, we mention that the Usadel equa- 
tions correspon ding to Eqs. (|6I7I) write, inside conductor 
Q (see e.g. Ref. Hfl) 

HDq §~z = [- t£f 3 + Mz) + iE ex (z)az,G] 

(El) 

The gap matrix A has a structure in the Nambu-spin 
subspace only, i.e., with our conventions, 



Using the angular parametrization of section [D] and 
A(z)eK, Eq.JED) leads to: 





A(z) 


2 ^ =( « + «rE 


A(z) = 


A(z) 




A(z)* 






A(z)* 





is + iaE ex (z)) sin(Ao-) — A(z) cos(A,j) 

(E2) 
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